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MAXIMUM LIFETIME ROUTING IN WIRELESS AD-HOC NETWORKS 
Field of Invention 

The present invention relates to wireless ad-hoc networks and, more particularly, 
to maximum lifetime routing techniques for use in wireless ad-hoc networks. 

5 Background of the Invention 

A wireless ad-hoc network, also referred to as a mobile ad-hoc network 
(MANET), is known to comprise a set of nodes connected by wireless links. Typical 
examples of ad-hoc networks are wireless sensor networks, where the nodes are sensors 
that gather environmental data and send the information to computational nodes for 

10 further processing, or to base stations for relay to a wired network. Such networks may 
be deployed, for example, in hazardous locations such as in disaster areas to aid rescue 
efforts, in areas for mineral or oil prospecting, and in battlefields for defense applications. 

The topology of an ad-hoc network is dynamic since nodes are typically free to 
move randomly and organize themselves arbitrarily. Therefore, the topology may be 

15 determined by the current geographic location of the nodes and other environmental 
conditions, and the characteristics of the radio transceivers that the nodes, possess. The 
topology may therefore be represented as an arbitrary graph with c< nodes" of the graph 
representing nodes in the network and "edges" of the graph representing links between 
nodes. 

20 The nodes in an ad-hoc network typically attempt to communicate amongst each 

other by relaying packets. However, due to the limited transmission range that is 
characteristic of nodes in an ad-hoc network, multiple network "hops" are typically 
needed for one node to exchange data with another node across the network. The 
problem is to design effective routing protocols to meet a variety of performance 

25 objectives given such a communications environment. 

Most existing routing protocols for wireless ad-hoc networks concentrate on 
finding and maintaining routes in the face of changing topology caused by mobility or 
other environmental changes. Typical protocols use shortest path methodologies based 
on hop count, geographic distance, or transmission power. The first two are important in 

30 minimizing delay and maximizing throughput. The third objective is peculiar to wireless 
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ad-hoc networks, and is important because typically the nodes involved have a limited 
power supply, and radio communication consumes a large fraction of this supply. 

To address this issue, several power-aware routing and topology control 
methodologies have been developed. In most of these methodologies, the aim is to 
5 minimize the energy consumed per packet in order to deliver it to the destination. The 
typical approach is to use a distributed shortest path methodology in which the edge costs 
are related to the power required to transmit a packet between the two nodes involved. 
The problem with this technique is that nodes on the minimum-energy path are quickly 
drained of power, affecting the network connectivity when they fail. 
10 Some of the more sophisticated routing methodologies associate a cost with 

routing through a node with low power reserves. But this remains, at best, a heuristic 
solution. 

A formulation using linear programming has been proposed which attempts to 
capture the issue of power consumption more precisely, see J.H. Chang et al., "Routing 

15 for Maximum System Lifetime in Wireless Ad-hoc Networks," Proceedings of 37th 
Annual Allerton Conference on Communication, Control and Computing, September 
1999; and J.H. Chang et al., "Energy Conserving Routing in Wireless Ad-hoc Networks" 
Proceedings of IEEE INFOCOM, pp. 22-31, March 2000, the disclosures of which are 
incorporated by reference herein. The idea is to make the goal of routing the 

20 maximization of the "network lifetime." As is known, the network lifetime is the time 
period in which the network operates without a node failure. This maximization approach 
utilizes a heuristic methodology to solve the linear program approximately. However, 
such a heuristic approach can perform arbitrarily poorly in worst case situations. 

Further, a centralized methodology to determine the maximum lifetime has been 

25 proposed in J.H. Chang et al., "Fast Approximation Algorithms for Maximum Lifetime 
Routing in Wireless Ad-hoc Networks," Lecture Notes in Computer Science: Networking 
2000, vol. 1815, pp. 702-713, May 2000, the disclosure of which is incorporated by 
reference herein. Such centralized methodology is based on the Garg-Koenemann 
methodology for multicommodity flow, as described in N. Garg et al., "Faster and 

30 Simpler Algorithm for Multicommodity Flow and Other Fractional Packing Problems, 
Proceedings of 39th Annual Symposium on Foundations of Computer Science, pp. 
300-309, November 1998, the disclosure of which is incorporated by reference herein. 
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However, this approach does not solve all of the performance deficiencies associated with 
the above linear programming approach. 

Accordingly, a need still exists for effective routing techniques that meet 
performance objectives associated with an ad-hoc network environment and the like. 

5 

Summary of the Invention 

The present invention provides effective routing techniques that meet performance 
objectives associated with an ad-hoc network environment and the like. The techniques 
of invention serve to substantially maximize the lifetime of the network. 

10 In one aspect of the invention, a packet routing technique for use in a node of a 

distributed network comprises the following steps/operations. Queues for storing packets 
are maintained, wherein at least one queue is associated with a link existing between the 
node and a neighboring node, and a queue has a height associated therewith. A route is 
then determined for one or more packets stored in the queues based on heights of queues 

15 at neighboring nodes, such that energy constraints associated with the node and the 
neighboring nodes are substantially maximized. As mentioned, the distributed network is 
preferably a mobile ad-hoc network wherein the node and at least one neighboring node 
communicate over a wireless link. 

One or more packets may be sent, over a link, from a queue of the node to a queue 

20 of a neighboring node when the height of the queue of the node is greater than the height 
of the queue of the neighboring node. Further, one or more packets may be received, 
over a link, in a queue of the node from a queue of a neighboring node when the height of 
the queue of the node is less than the height of the queue of the neighboring node. 

The node may receive broadcast information from at least one neighboring node 

25 pertaining to the height of at least one queue of at least one neighboring node. Further, 
the node may broadcast information to at least one neighboring node pertaining to the 
height of at least one queue of the node. 

Packets may be routed to and/or from the node in rounds such that throughput 
requirements are substantially satisfied while substantially maximizing a lifetime 

30 associated with the distributed network. Further, the network lifetime has an upper bound 
and a lower bound associated therewith. 

Further, route determination techniques of the invention may also take into 
account one or more of idle power consumption associated with a node, computation 
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power consumption associated with a node, a periodic recharge associated with a node, 
one or more edge constraints, and power consumption associated with receiving a packet 
at a node. 

Advantageously, as will be explained in detail herein, the invention may provide a 
5 distributed routing methodology that reaches an optimal routing solution to within an 
asymptotically small relative error. Further, the techniques of the invention are preferably 
based on the formulation of multicommodity flow, and they take into consideration 
different power consumption models and bandwidth constraints. While, the inventive 
techniques may be employed with a variety of networks, they work particularly well for 
10 both static and slowly changing dynamic networks. However, the inventive techniques 
may also be implemented in accordance with a dynamic network. 

These and other objects, features and advantages of the present invention will 
become apparent from the following detailed description of illustrative embodiments 
thereof, which is to be read in connection with the accompanying drawings. 

Brief Description of the Drawings 

FIG. 1 is a block diagram illustrating a wireless ad-hoc network in which the 
techniques of the present invention may be implemented; 

FIG. 2 is a block diagram illustrating an example of a computing system which 
may be employed to implement the techniques of the present invention; 

FIG. 3 is a flow diagram illustrating a core routing methodology according to an 
embodiment of the present invention; 

FIG. 4 is a flow diagram illustrating a distributed routing methodology according 
to an embodiment of the present invention; and 

FIG. 5 is a flow diagram illustrating a decision flow at each node for a distributed 
routing methodology according to an embodiment of the present invention. 

Detailed Description of Preferred Embodiments of the Invention 

The following description will illustrate the invention using an exemplary ad-hoc 
network architecture. It should be understood, however, that the invention is not limited 
to use with any particular ad-hoc network architecture. Rather, the invention is more 
30 generally applicable to any network architecture in which it is desirable to perform 
effective data routing. 
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Referring initially to FIG. 1, a block diagram illustrates a wireless ad-hoc network 
in which the techniques of the present invention may be implemented. As shown, a 
wireless ad-hoc network 100 comprises a plurality of mobile nodes 102-1 through 102-N 
connected by wireless links. The number of nodes N depends on the application in which 
they are to be deployed and, thus, the invention is not limited to any particular number. 

Nodes 102-1 through 102-N are free to move randomly and organize themselves 
arbitrarily. Thus, the topology of the ad-hoc network 100 is dynamic. Further, nodes 
102-1 through 102-N communicate amongst each other by relaying packets. However, as 
mentioned above, due to the limited transmission range that is characteristic of such 
nodes in an ad-hoc network, multiple network hops (denoted as 104 in FIG. 1) may be 
needed for one node to exchange data with another node across the network 100. 

Thus, as shown in FIG. 1, node 102-1 may communicate with node 102-N by 
transferring packets via a wireless link to node 102-3 (first hop), which then transfers the 
packets via a wireless link to node 102-4 (second hop), which then transfers the packets 
via a wireless link to node 102-N (third hop). Of course, depending on various factors, 
packets may be transferred from node 102-1 to 102-N via other routes with more or less 
hops 104. It is also to be understood that, for the sake of clarity, FIG. 1 does not 
expressly illustrate all possible hop arrangements that may be utilized to transfer packets 
from one node to another. Also, it is to be appreciated that a node may have separate 
wireless links between itself and other nodes, however, a hop refers to the transfer of data 
from the node to one of the other nodes via one of its wireless links. In any case, it is to 
be appreciated that it is the selection of a route (i.e., determining a routing solution) in the 
network upon which the invention primarily focuses. 

It is to be appreciated that any one of the nodes in the network 100 may serve as a 
computational node, wherein the packets transferred thereto may be further processed in 
some manner. In addition, node 102-N, or any other receiving node within range of a 
base station 106 in the network, may transfer information to the base station for relay to a 
wired network 108. The wired network 108 may be a private data network, the Internet, 
etc. It is to be understood that, for the sake of clarity, only one base station and wired 
network is illustrated in FIG. 1; however, there may be one or more additional base 
stations communicating with one or more additional wired networks. Of course, it is also 
to be understood that network 108 could alternatively be another wireless network. 
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Thus, in the case where the nodes 102-1 through 102-N are deployed in a 
hazardous or hostile environment in order to collect data in the environment, the data 
collected by the nodes may be transferred to a remote location, via the wired network 108, 
such that it can be presented and/or analyzed in a non-hazardous or non-hostile (or, at 
5 least, less hazardous or hostile) environment. 

Referring now to FIG. 2, a block diagram illustrates an example of a computing 
system or device which may be employed to implement all or a portion of the techniques 
of the present invention. More particularly, one or more of nodes 102-1 through 102-N 
and/or base station 106 may implement such a computing system 200 to perform the 

10 techniques of the invention (although, a base station may not necessarily require a 
"sensor," as will be explained below). Of course, it is to be understood that the invention 
is not limited to any particular computing system implementation. 

In this illustrative implementation, a processor 202 for implementing at least a 
portion of the methodologies of the invention is operatively coupled to a memory 204, a 

15 sensor 206 and a network interface 208 via a bus 210, or an alternative connection 
arrangement. It is to be appreciated that the term "processor" as used herein is intended 
to include any processing device, such as, for example, one that includes a central 
processing unit (CPU) and/or other processing circuitry (e.g., digital signal processor 
(DSP), microprocessor, etc.). Additionally, it is to be understood that the term 

20 "processor" may refer to more than one processing device, and that various elements 
associated with a processing device may be shared by other processing devices. 

The term "memory" as used herein is intended to include memory and other 
computer-readable media associated with a processor or CPU, such as, for example, 
random access memory (RAM), read only memory (ROM), fixed storage media (e.g., 

25 hard drive), removable storage media (e.g., diskette), flash memory, etc. The memory 
may preferably be used to implement the one or more queues associated with a node, to 
be discussed below. 

Further, the term "sensor" as used herein is intended to include, for example, one 
or more devices capable of collecting data from the environment in which the computing 
30 system 200 (e.g., mobile node) is deployed. Examples of sensors may be receivers 
associated with wireless communication devices, which are capable of collecting data 
from the environment in which they are deployed. However, it is to be understood that 
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the invention is not limited to any particular sensor. Also, a base station will typically not 
have a sensor. 

Still further, the phrase "network interface" as used herein is intended to include, 
for example, one or more devices capable of allowing the computing system 200 (e.g., 
5 mobile node or base station) to communicate with another computing system (e.g., 
mobile node or base station). Thus, the network interface may comprise a transceiver 
configured to communicate with a transceiver of another computing system (e.g., mobile 
node or base station) via a suitable wireless protocol. Examples of wireless transmission 
protocols may be Bluetooth and IEEE 801.11. However, it is to be understood that the 

1 0 invention is not limited to any particular communications protocol. 

It is to be appreciated that while the present invention has been described herein in 
the context of a network communications system, the methodologies of the present 
invention may be capable of being distributed in the form of computer readable media, 
and that the present invention may be implemented, and its advantages realized, 

15 regardless of the particular type of signal-bearing media actually used for distribution. 
The term "computer readable media" as used herein is intended to include recordable-type 
media, such as, for example, a floppy disk, a hard disk drive, RAM, compact disk (CD) 
ROM, etc., and transmission-type media, such as digital and analog communication links, 
wired or wireless communication links using transmission forms, such as, for example, 

20 radio frequency and optical transmissions, etc. The computer readable media may take 
the form of coded formats that are decoded for use in a particular data processing system. 

Accordingly, one or more computer programs, or software components thereof, 
including instructions or code for performing the methodologies of the invention, as 
described herein, may be stored in one or more of the associated storage media (e.g., 

25 ROM, fixed or removable storage) and, when ready to be utilized, loaded in whole or in 
part (e.g., into RAM) and executed by the processor 202. 

In any case, it is to be appreciated that the techniques of the invention, described 
herein and shown in the appended figures, may be implemented in various forms of 
hardware, software, or combinations thereof, e.g., one or more operatively programmed 

30 general purpose digital computers with associated memory, implementation-specific 
integrated circuit(s), functional circuitry, etc. Given the techniques of the invention 
provided herein, one of ordinary skill in the art will be able to contemplate other 
implementations of the techniques of the invention. 
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For ease of reference, the remainder of the detailed description is divided into the 
following sections: (1) Introduction; (2) Problem Formulation; (3) Local-Control Flow 
Methodology; (4) Maximi zi ng Lifetime; (5) Implementation; (6) Extension to Local 
Constraints; and (7) Illustrative Decision Flow at Node. 

5 1. Introduction 

As mentioned above, given the characteristics of an ad-hoc network environment, 
a problem is to design effective routing methodologies to meet a variety of performance 
objectives of the ad-hoc network environment. The invention addresses this problem and 
provides solutions for realizing such performance objectives. 

10 In accordance with the invention, as illustrated hereinbelow, the routing problem 

is explicitly formulated as a maximum concurrent flow problem. All routing problems 
have analogous flow realizations, where the flow represents the routes that packets take, 
and the demands represent the rate at which packets are generated by various nodes. This 
formulation allows for the application of maxflow techniques to the problem of routing 

15 for maximum network lifetime. However, many of the known maxflow techniques are 
not well-suited to distributed implementation as is necessary for a routing methodology, 
or do not adapt well to changing network topology. Moreover, they are typically 
better-suited to the problem of wired network routing where the constraints are link 
bandwidth (the number of packets that can be routed over an edge in unit time) and node 

20 capacity (the number of packets that can be forwarded by a node in unit time). 

The present invention provides a distributed flow methodology to the routing 
problem. The methodology finds an approximation to a feasible flow if one exists. 
Advantageously, unlike existing heuristic approaches (e.g., heuristic approaches described 
above in the background section), the approximation factor of the invention is guaranteed. 

25 Also, in accordance with the invention, lower bounds are established for the performance 
of the methodology. While the techniques of the invention are applicable to a wide 
variety of networks, the techniques are particularly well-suited to static networks with 
various power consumption models, as well as dynamic networks with slowly changing 
edge costs. 

30 One example of an advantage that the inventive approach has over existing routing 

approaches is that the invention is a distributed, local-control approach and therefore does 
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not require a central node with global knowledge of the network. This is also one reason 
why the inventive approach works well with dynamic networks. 

2. Problem Formulation 

In this section, the problem is formulated more precisely, and both a linear 
5 programming and a flow interpretation are provided. 

The subject network in which the inventive routing methodology is to be 
implemented (such as, for example, network 100 in FIG. 1) is represented as a graph, with 
N nodes and M edges, with the nodes representing wireless devices and the edges the 
wireless links between them. Associated with each node i is a quantity E i9 representing 
10 the initial energy reserve of the device. Each edge ij has a cost which is the energy 
required to transmit one packet of data across the corresponding link. 

The routing problem comprises of a set of K source and destination pairs, which 
are pairs of nodes in the network. Each pair has a throughput requirement Q Cy the number 
of packets per second that must be routed between source and destination nodes. The 
15 routes are given by variables^, which represent the rate at which packets are transmitted 
across link ij for connection c. 

2.1 Linear Programming Formulation 

The "lifetime" of the network is defined as the time period before which the first 
node failure occurs, that is, the time period before which some node's energy reserve is 
20 reduced to zero. The invention routes packets in such a maimer that the lifetime is 
maximized while the throughput requirements are satisfied. The lifetime is denoted by T. 
Variables/^ denote the total number of packets for connection c transmitted from node i 
to node j over the lifetime of the network. 

The total energy consumed at node i is given by: 

25 2 ea # 

where the sum is over all nodes adjacent to i and all connections c. At every node except 
the source and destination for a particular connection c, the number of packets received 

9 
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equals the number of packets transmitted. At the source, the number of packets 
transmitted equals Q C T over the lifetime T of the network. 
Hence, the linear program is: 

maximize T 

such that: (1) fg^O Vy',c ' 
(2) V/ 




i is a source for c otherwise. 



The second constraint is due to the finite power supplies at the nodes, and the third 
represents the throughput requirements at the sources and the conservation constraints at 
other nodes. 

This formulation serves to show that the problem can theoretically be solved 
exactly in polynomial time. 

2.2 Multicommodity Flow 

The routing problem can generally be interpreted as a flow problem. The 
variables can be thought of as flow values of a commodity c. The requirement that 
incoming packets arrive at the same rate as outgoing packets leave is equivalent to the 
flow conservation constraint that, at any intermediate node, the amount of incoming flow 
equals the amount of outgoing flow. The sources and destinations of the routing problem 
become sources and sinks in the flow problem. The throughput requirements Q c represent 
the rate at which flow is produced at the sources and consumed at the sinks. So at each 
instant of time, the routing is represented as a multicommodity flow satisfying demands 
Qc 

The time evolution of the network is considered as proceeding in a sequence of 
synchronous rounds, each lasting for unit time. In each round, Q c packets are generated at 
each source node. Packets may be transmitted between nodes, and the energy reserve of 
each node is reduced by the total cost of the packets it transmits. So, iffi is the number of 
packets transmitted by node i to node j 9 the energy reserve of i is reduced by Sy eyfy in 
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this round. Since each node will actually be forwarding packets for a number of 
connections, the reduction in the energy is given by: 

From the problem formulation, it might appear that a more general version is to 
5 ask for how many rounds a feasible multicommodity flow can be maintained, where the 
flow is allowed to vary with time, and some intermediate nodes might fail before the 
demands become infeasible. This does not lead to an improvement in the lifetime 
achievable over a static flow that lasts only until the first node failure, however, as can be 
seen by the following argument: suppose the flow achieving the maximum lifetime T is 
10 fij CO at time U and satisfies the demands at every instant. Replace this variable flow by 
its time average: 

j>= \llfft)dt 

This static flow will satisfy the same demands at every instant, and now no node 
will fail before the lifetime T. Since the flow is static, the energy consumed in any round 
15 at a node i cannot exceed E ilT. 

Thus, it is enough to look for a multiflow satisfying all the demands, and 
constrained by: 

2e^/|^ 4r for all i 

to achieve a lifetime of T. An approximation methodology for this problem will be 
20 described in the next section. 

3. Local-Control Flow Methodology 

The methodology of the invention will be expressed in terms of the 
above-described flow interpretation. A goal is to determine whether a given value of the 
lifetime T is feasible. That is, is it possible to route flow in the network so that all the 
25 demands are satisfied, and no node runs out of energy before T rounds have gone by. 

11 
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It is assumed that there exists a multiflow that satisfies not just the demands Q C9 
but slightly more, namely, the demands (l+e)Q C9 for some parameter e > 0. 

Thus, the present invention provides the following principle (1): 

If the demands (l+e)Q c are feasible, then there exists a routing methodology that 
5 satisfies the demands (1- 5)Q C , for any 8>0, after running for ^y Wg)) rounds. 

Principle (1) will be proven below by actually describing such a methodology and 
giving a complete analysis of its performance. It is to be understood that, in accordance 
with the inventive approach, the constraints are on the nodes, since each node must 
respect its energy budget E i/T. 

10 3.1 Description 

Each node maintains K queues for each of its links. In each round, the flow in a 
queue may be moved to the other end of the link. Flow may also be redistributed among 
the queues at a single node. 

Denote by qjj the height of the queue for commodity c on link ij at node i. The 
1 5 link has another queue at its other end, denoted by qf t . 

A potential function <f>(q) associated with every queue is defined, where q is the 
queue height. The potential function will be completely specified when the performance 
of the methodology is analyzed below, for now note that it is a twice-differentiable 
convex function. The methodology operates so as to minimize the sum of the potential 
20 functions of all queues. 

Referring now to FIG. 3, a flow diagram illustrates a routing methodology 
(referred to herein as the "core" methodology) according to an embodiment of the 
invention. In each round of the methodology 300, the following phases (steps) are 
performed: 

25 • Inject flow (step 302). For each commodity c, add Q e flow at the corresponding 
sources. 

• Balance nodes (step 304). Equalize the queues at each node. 

• Push flow (step 306). Associate each pair of queues q§ and qf t for an edge ij with the 
node where the queue height is greater. To reduce the potential, packets are moved 

30 from the higher queue to the lower (this follows from the convexity of 0). So, the 

following function is minimized: 

12 
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where the sum is over edges if associated with node i 9 and the variables are the/£ subject 
to the constraint: 

The solution is then used to move flow from node i to nodey . 

• Drain flow (step 308). Absorb the flow for commodity c at its sink. That is, set the 
height of all queues for c that are at the sink to zero. 

3.2 Analysis 

It will be shown that the potential in the system remains bounded as time goes on, 
which implies that most of the flow injected into the network reaches the appropriate sink. 

The idea is to lower-bound the decrease in potential in a round by comparing what 
the methodology does to what would happen if the flow values were picked according to 
a feasible multiflow instead. This will in turn imply an upper bound on how large the 
potential can become as the methodology is repeatedly executed. 

For simplicity in writing the expressions, it will be assumed that the source nodes 
and the destination nodes both have unit degree, and that all flows and queues are scaled 
so that the demand for each commodity is unity. 

The increase in potential after the first two phases will be at most: 

2^(5 C )-^(5 C -1)^E^ , (5 C ) ' 

where s c is the total amount of commodity c at the source after flow injection. 

The potential decreases when flow is pushed, and to estimate the decrease in 
potential, it is considered what happens if a feasible flow// satisfying the demands (1 + e 
) is used, instead of the computed minimizing flow values. The expression is: 
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<sr A) + * to) - <t> to +/20 - 4> to -m 

where the summation extends over all edges in the network, and q% and q£ denote the 
queue heights for commodity c at the head and tail of the edge, respectively. The 
extended mean value theorem is applied to <j> and the expression simplified to: 

- iCTS) 2 + AWT) + r (g y e e )) 

where the Ce/ and are suitable constants between 0 and 1. Note now that the first 
derivative terms will cancel out except at the source and sink, because in the node 
balancing phase, all the queue heights for a given commodity at a given node are made 
equal, and// satisfies the flow conservation property. The potential function is chosen so 
that the second derivative terms make a very small contribution to the sum. In particular, 
assume that their contribution is at most f <j> ' ( Sc ). Hence, the potential drop is at least: 

L(l+e)((l- fty'tec)- *>'((») 
So, overall during the round, the potential drops by at least: 
?(f" f>0c)-(l + e)^(0) 

Note that if the flow values ff that move flow from a lower height queue to a 
higher queue are ignored, only the potential drop will be improved. Since, the inventive 
methodology considers all such feasible updates while minimizing, the potential drop it 
finds will be better than the one just derived. 

Let S be the value that satisfies the equation: 

(f- f>'CS)= (l+eW(0) 
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So, if s c £ S for any commodity, then the potential drop during the round will be 
non-negative. The amount of commodity c in the network is bound by IMS by showing 
that no queue can exceed S in height. 

To accomplish this, overflow buffers are used at each source node. The idea is 
that if the queue at the source has reached S, then additional injected flow is put in a 
special overflow buffer, which has potential b<f> ' (5) if it has b amount of commodity. 
With this potential function, the increase in potential when flow is injected remains 
upper-bounded by <j>\s c ). The advantage of bounding the source queue is that no other 
queue can exceed this bound. It is thus easy to see that, given a particular value of queue 
height, S in this case, the first queue to exceed this height must be at the source. 

Now we argue by induction on the number of rounds that the potential must 
always be less than 2MK<f> (S). For as long as the overflow buffers are empty for all 
commodities, this remains true. If, on the other hand, the flow injection phase overflows 
some source queues, then the potential actually drops during this round. 

Hence, the maximum size of an overflow buffer is 2MK<j> (S)l<f> ' (S). This result 
will be needed in the next section, so it is restated as principle (2): 

IfS satisfies the equation 

( f - -f>'CS)= (l+e)K<pi(0) 
then the maximum size of an overflow buffer is 

2MK<(> WCS). 

The maximum amount of commodity c within the edge queues is IMS, hence, 
combining this with the upper bound on the overflow buffers, we get an upper bound of 
2M{S+K<f> (S)/ft(S)) on the total amount of commodity in the system. 

Now we set the potential function </>{q) = e^,foxa= e/AL, where L is an upper 
bound on the length of the longest flow path in the feasible flow//. This gives, for e 
small enough: 

? We e ) 2 (.t»(q e eH + C%fe C ) + <K<& " CW» 
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^(l + e)ae a(1+€) 2^ / (^) 
<; (1 + s)^e ■ ii^V4L L(1 + 

5 . The equation for 5 gives ^( w * We) ). So, the total amount of a given 

commodity in the system is at most R = 8ML(K+ ln(6K/e))/e. 

Suppose that the methodology is run continuously for RJd rounds. The amount of 
commodity c pumped into the network is R/S 9 and the amount remaining in the network is 
at most R. Hence, the average amount of commodity transported through the network per 
10 round is 1 - 8. This proves principle (1) above. 

Another parameter of interest is the average delay a packet experiences while 
being routed by the inventive methodology. Since the amount of commodity in the 
system is bounded by R in the long run, the average delay is also upper-bounded by R 9 by 
Little's Law. 

1 5 Thus, principle (3) may be stated as: 

The average delay of a packet is at most ^K+uiWe)) 

4. Maximizing Lifetime 

Now a description is provided on how to use the methodology to converge to the 
routing appropriate to the maximum lifetime. Consider temporarily that centralized 
20 control is available, and information can be updated at all the nodes at the same time. 

The node constraints will be set to Ei/T, and we will vary T until it converges to 
the optimal value T. The idea is that if T is too large, then the node constraints are too 
tight to allow for a feasible flow, hence the source queues will grow without bound. 
Hence, we start out T at some value that is smaller than the true lifetime and increase it 
25 until we find that some overflow buffer exceeds the bound 2MK<f>($)l<f> \S). Now, we 
know the optimal lifetime T to within a factor of two, and can perform a bisection search 
to determine it exactly. 

Suppose it is desired to determine the optimum lifetime T to within a relative 
error of e. Initialize T to some suitable value r 0 , which is assumed to be less than the 
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optimum T. Let the methodology run for 2R rounds. If no overflow buffer exceeds the 
bound 2MK<f>($)/<f> f (S\ then it is known that the average fraction of commodity 
transported through the network is at least one-half by principle (1). 

Notice that if demands XQ C can be satisfied with constraints E t /T 9 then scaling all 
the flow values by XIX shows that demands Q c can be satisfied with constraints E t IXT. 
Hence, if demands QJ2 are being satisfied on average with constraints E { /T, then it is 
known that demands Q c can be satisfied with constraints 2E ( /T 9 and so T* ^ 772. T is 
doubled and run for another 2R rounds. This process is repeated until some overflow 
buffer bumps up against the 2MK<f>(S)l<f> '(£) upper bound, say when r= T x . At this point, 
it is known that: 

(i+e)ri^r*^^- 

The upper bound is because if demands (l+e)<2c had been feasible with 
constraints E t ITu then by principle (2), the overflow buffers cannot exceed 
2MK<f>(S)/<f>'(S). 

The methodology is stopped as soon as an overflow buffer hits its bound, even if 
the full 2R rounds have not been executed. Because of this, the maximum overshoot that 
can happen in an individual buffer is one unit, corresponding to a potential overshoot of 
<f> '(5). To get the potential back under the bound 2MK(j>{S\ a round of the methodology is 
run without injecting any new packets at the sources that overshot. Suppose that K 1 of 
the source overflow buffers exceed their bounds. Then, the decrease in potential that this 
round produces will be at least K'j'QS), and hence will be sufficient to reduce the 
potential below 2MK<f>(S). 

The power consumed at node i during this process of finding an upper bound on 
T* is at most: 

Now that upper and lower bounds are obtained, a search is performed to determine 
the value of J* more accurately. The first approach is to use a bisection search, running 
the methodology for RIs rounds to cut the search interval approximately in half. But a 
second technique cuts the search time significantly. If T h and T\ are the current upper and 
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lower bounds (initially T h = Tx and 77 = 7V4), a test value T= (277> + 7/)/3 is chosen, a 
third of the way from T h to77, and the methodology is run for Rid rounds, where 
<5 = (Za - Ti)/3T. The idea is that if no source overflow buffer exceeds its bound within 
these many rounds, then by principle (1): 

5 (1 + e)T h >T*l>7{l-d) = 

and if some buffer does exceed its bound, then: 

(l+e)T^T*>Ti 

In either case, the search interval is reduced by a factor of 2/3. The excess power 
consumed at node i (over the ideal E t /T*) during these Rid rounds is at most: 

10 KC^l Il\ ^ 3E i R 

The number of rounds spent during this search process is at most: 



and the total excess power consumption at node i is at most: 

15 3£//glog 3/2 (3/e) 

The total number of rounds spent to converge to an accurate value of T* is hence: 

<21og4jP + #) 

During the remaining rounds, the power consumed per round is at most 
20 Eil{\- e)T*. Hence, if: 
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AEjR 3£f/?log^(3/g) 
To + 7^ ^ 

then the node will last for at least (1 - 2e)T* rounds. Hence, if To = SR/e is chosen, and 
if: 

5 then the remaining life of the node is at least (1 - 2 e )r*. 
Further, a principle (4) may be stated: 

There exists a routing protocol that achieves a lifetime of (I - 2e)T*, if 
For a given value of T*: 

io G = ©(^ + io g ^))i 

should be chosen to get the best possible bound. It is to be understood that /= ©(g) 
means that there exist constants c\ and cz such that c\g <.f<. czg. 

5. Implementation 

In this section, three issues associated with implementing the methodology as 
15 presented so far are considered. The first issue is that, thus far, the flow has been treated 
as a continuous quantity. However, in routing, it is typically not desirable to split up 
packets. Thus, the queues are discretized so that they contain an integral number of 
packets, and an entire packet is routed at a time. The second issue is computing the flow 
values that minimize the potential function at a node. The third issue is how to converge 
20 to the maximum lifetime in a distributed manner rather than in the centralized fashion 
analyzed above. 

Another problem that occurs in practice is that traffic does not arrive in a smooth, 
deterministic maimer, but may be intermittent or bursty in nature. To deal with this, a 
"leaky bucket* * is added at every source to smooth out the traffic. 
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5.1 DiscretMng Flow 

The node rebalancing phase of the methodology is modified to rebalance the 
queues while maintaining an integer number of packets in each queue. This implies that 
the actual height of a queue might differ by up to one packet from its ideal height. 

This introduces an error term into the analysis of the methodology. While 
deriving the drop in potential by routing a flow f* across an edge e, the following 
expression is obtained: 

2 fSdfi '(q%t) - <t> '(flfo)) + second derivative terms 

The first term might differ from its ideal value (when the q's are allowed to vary 
continuously) by: 

£ fM "fa + tfJQ.) + 4> "(gtH + 

where 77^ and tj c eh are between -1 and 1. This sum is at most: 

g §(2a 2 e <*) 2(1 + e)L^ '(S) <; e<f> '(S) 

So, this does not significantly affect the analysis of the methodology. 

5.2 Minimizing Potential 

The problem to solve at each node is: 

Minimize }jin C ii-fv) + Klji+fi) 

subject to 

Using the technique of Lagrange multipliers, the/^ satisfy the system: 
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So, the maximum of the quantities is computed: 

e u 

5 (which must be equal when the minimum potential solution is reached) over the edges ij 
and connections c and route one packet for that edge and connection. The queue heights 
are updated and the process repeated until the energy budget EtlT 'is exhausted. 

5.3 Distributed Implementation 

The core of the methodology (described above in section 3) may be implemented 
10 in a distributed manner in a straightforward manner. A more difficult task is defining a 
distributed protocol for doing the work of maximizing the lifetime (described above in 
section 4). Such a distributed protocol is described below. 

The idea is that each source sends out a broadcast (forwarded to each neighbor 
along with data packets) whenever its overflow buffer hits the bound 2MK<f>(S)l^ ! (S). 
1 5 Each node receives the broadcast in at most N rounds. 

In the initial phase, every 2R rounds, every node doubles its local value of T. If 
during the next N rounds, it receives a broadcast message from a source, it switches into 
the next phase, and forwards the broadcast. 

In the second phase, each node maintains the upper and lower bounds Th and 77, 
20 and computes the test value T= (27/, + 77)/3 to determine its local energy budget. If after 
RI5 + Grounds, no source broadcast is received, the value of 77 is updated to (7* + 277 )/3. 
If a broadcast is received, the value of Th is updated to T This continues until the value 
of T* is determined (to within relative error e). 

Referring now to FIG. 4, a flow diagram illustrates a routing methodology 
25 (referred to herein as the distributed routing methodology) according to an embodiment of 
the invention. The methodology 400 executed at each node can be summarized as 
follows: 

• Initialize (step 402). Setr=7o. 
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• Determine upper bound Th (step 404). Run core methodology until 2R + N 
rounds pass, a broadcast is received, or (if the node is a source) the overflow 
buffer hits 2MK<f>(S)/<f>'(S). If 2R + N rounds have gone by, double T and 
repeat this step, otherwise move to the next one, 

5 • Broadcast (step 406). If a broadcast message was received, forward it to 

neighbors, else generate a new one. 

• Search initialization (step 408). Set T h = Tand 77 = 774. 

• (step 410). Compute T= (27* + 77)/3,<5 = (T h - 77)/3. Run the core 
methodology until Rid + Grounds pass, a broadcast is received, or (if the node 

10 is a source) the overflow buffer hits 2MK^(S)/(ff'(S). 

• Update (step 412). If R/3+N rounds have passed, set 77 = (7* + 277)/3. 
Otherwise, setT* = (27* + 77)/3 and either forward the broadcast or generate a 
new one. If {T h - 77)/77 £ e, go back to the Search step. 

The analysis of this distributed protocol is similar to the centralized one and leads 
15 to the same bounds except for a small constant factor, since N is small compared to R. 

5.4 Arbitrary Traffic Sources 

In practice, traffic may not arrive in the smooth, deterministic maimer we have 
assumed so far, but may be intermittent or bursty in nature. To deal with this, we add a 
leaky bucket at every source to smooth out the traffic. A leaky bucket is a controller used 
in network applications to smooth out traffic flow. The idea is that p tokens are generated 
per second, and tokens are dropped if there are more than a of them. Each packet 
transmission consumes one token. If tokens are not available, the packets are buffered 
until they become available. Such a controller limits the traffic transmitted: the number 
of packets transmitted in a time interval T can never exceed pT+ a. In the present 
application, the special case where a is zero is employed, so that the traffic is smoothed 
out completely. 

5.5 Acceleration of the Methodology 

It is realized that a delay can grow unexpectedly fast. On a route of length L, the 
delay a packet faces may grow as L 2 because the queues decrease linearly from source to 
30 sink (after convergence to the appropriate flow). To both reduce this delay and to 
accelerate convergence, the methodology may be modified in the following maimer. The 
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flow to be sent on an edge may be modified by using some history, that is, the actual flow 
is computed as a linear combination of the flow that the original methodology (or the 
first-order method) chose, and the flow that was sent on the edge in the previous step. For 
a suitable choice of the coefficients of the two components, the convergence can be 
accelerated in this new second-order method. If the coefficient of the flow in the previous 
step is close to 1, then the final queue heights will be smaller than in the first-order 
method, thus reducing delay. 

The actual flow in the second-order method is thus calculated as: 

f=afo+pf 

where /o is the flow given by the first-order method,/' is the flow sent in the previous 
step, and a and ft are suitable parameters. The convergence can be significantly 
accelerated when a + > 1 . The closer $ is to 1 and the larger the value of a, the smaller 
the final queue heights and the faster the convergence. However, if a is too large, the 
methodology may become unstable. 

An exact analysis can be done for the simple case of a path graph, which shows 
that the second-order method is stable as long as a < 1 + /?. 

This method may call for flow values that are greater than the actual queue 
heights. Two methods to deal with this include (i) sending the maximum possible when 
this happens, and ignoring the rest; or (2) the excess flow called for can be remembered 
and sent during subsequent time steps if possible. 

5.6 Alternative Implementation and Simulation Results 

In a simplified, practical implementation of the methodology, the potential 
function may be taken to be just the square of the queue heights rather than the 
exponential function. 

Also, instead of using a distributed protocol to search for the true lifetime, the 
potential function may be modified by adding a term proportional to the square of the cost 
of the flow, i.e.: 
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This works very well as long as the constant a is chosen appropriately. With this 
modified potential function, choosing the flow values to minimize the potential becomes 
reasonably simple to do quickly and exactly. 

For the simulation, the graphs were generated by randomly generating points in a 
5 square area. Performance was measured by first computing the maximum lifetime using 
the methodology from J.H. Chang et al., "Fast Approximation Algorithms for Maximum 
Lifetime Routing in Wireless Ad-hoc Networks," Lecture Notes in Computer Science: 
Networking 2000, vol. 1815, pp. 702-713, May 2000, running both the first- and the 
second-order methods for the hfetime number of rounds and noting how much flow 
10 actually reached the sink. 

The simulation shows that choosing a appropriately is largely a matter of making 
sure that it is not too small. Larger values reduce the total flow delivered, but only 
slightly. 

6. Extension to Local Constraints 

The inventive methodology works for a wide variety of different constraints, 
either on the nodes or on the edges or even both. 

From the proof of the core methodology given above, it follows that the 
methodology works for any combination of local constraints (both nodal as well as edge 
capacity constraints) as long as there exists a flow satisfying demands 1+e and meeting 
all the constraints. These constraints need not even be linear, as long as we are able to do 
the potential minimization. 

6.1 Idle Power 

Typically, communication is not the only energy consumer in a wireless node. 
The node may be doing some computation or data collection which consumes power at a 
rate independent of the number of packets it transmits. Also, wireless interfaces consume 
power even when idle and not receiving or transmitting any packets. To take these factors 
into account, a constant power consumption C t may be added to the node energy 
constraint, which becomes: 
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The core methodology is not affected at all by this change. The only thing that 
need be verified in order to show that the lifetime maximization protocol will still work is 
the scaling assumption that if demands XQ C can be satisfied with lifetime T, then demands 
Q c can be satisfied with lifetime XT. This follows because if flow values f§ achieve 
demands XQ C9 then the flow/^/A will satisfy Q C9 and: 

provided X <> 1, which is true whenever this result is applied {X = 1 -S in section 4). 

6.2 Periodic Recharge 

An interesting case is when some nodes are recharged periodically, say every Tt 
seconds. In this case, the right hand side of the constraint for this node can be replaced 
with E t ITi, since the node only has to survive until its next recharge. This node does not 
participate in the lifetime-determining phase of the protocol, except for forwarding 
control packets. 

6.3 Edge Constraints 

Edge constraints, for example, bandwidth constraints, can be incorporated. At 
each node, in addition to the energy constraint: 

additional constraints: 

may be included where By is the maximum possible rate at which packets can be 
transmitted on link ij, perhaps due to signal-to-interference ratio (SIR) requirements. 
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It may also be necessary to modify the energy constraint to allow non-linear 
power dependences. For example, there may be an initial startup cost associated with 
using the transmitter or for setting up the wireless link with the node at the other end, 
which may be a significant fraction of the total power required. This may be 
5 implemented in accordance with the inventive methodology. 

6.4 Receive Power 

A transceiver consumes power when operated in receive mode as well as in 
transmit mode. If the network is fairly static or slowly evolving, then this can be taken 
into account by charging it when the packet is actually forwarded instead of when it is 
10 received. Since all but a bounded number of packets at any node will be forwarded, the 
energy required to receive a packet is added to the energy required to transmit it in the 
node constraint, which takes the form: 

where e\j is the energy per packet for transmission, and e£ is the energy per packet for 
1 5 reception and decoding. 

6.5 Dynamic Network 

The core methodology described in section 3 works even if the network is 
dynamic and edge costs are varying with time, as long as it is true that there exists a 
multicommodity flow satisfying the demands Q c . 

20 However, because of the changing edge costs, it is possible that the lifetime may 

change. If any source overflow buffer hits its bound, then the current value of T is too 
high, and a new search can be begun using the distributed protocol with T h = T and 
Ti = 0. If the current value of T is too low, this will not be detected unless the entire 
lifetime maximization protocol is reinitiated. This can be done after at least R/e rounds 

25 have gone by, using for the starting value To the current value of T. Waiting for this 
many rounds ensures that 1 -e fraction of packets are transmitted to the appropriate 
destinations. So, if the rate of change in the network is slower than this, the lifetime may 
be tracked as the network evolves. 
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7. Illustrative Decision Flow at Node 

Accordingly, as illustratively explained above in sections 1 through 6, the 
invention analogizes network traffic to fluid flow. Using such an analogy, data buffers at 
nodes may be thought of as fluid vessels, wherein buffer height corresponds to pressure. 
As is known, fluid tries to equalize heights. Thus, in the context of network traffic, the 
techniques of the invention make routing decisions based on local information about 
queue heights at neighbors. 

Advantageously, at each node, queues are equalized over edges. Packets are 
routed from higher queues to lower queues, for each edge. Such routing is constrained by 
an energy budget, such that less data is routed over high-cost edges. Further, there is a 
trade-off between edges and connections, and such a trade-off is controlled by a potential 
function (see, e.g., <f> discussed above). 

Referring now to FIG. 5, a flow diagram illustrates a decision flow at each node 
for a distributed routing methodology according to an embodiment of the invention. The 
decision flow 500 starts at block 502. At step 504, the node determines if a packet has 
arrived. If yes, the node determines, in step 506, if the packet is for it. If yes, the packet 
is accepted, in step 508, and the process returns to step 504. If the packet is not for the 
node, the packet is put in a local queue, in step 510. The process returns to step 504. 

Returning to the initial determination in step 504, if no packet arrives, the node 
determines whether there are any packets in its local queue, in step 512. If no, the process 
returns to step 504. If yes, in step 514, the node determines if there is a neighboring node 
that has a shorter queue. If no, the process returns to step 504. If yes, the packet is sent to 
the neighbor, in step 516. 

Although illustrative embodiments of the present invention have been described 
herein with reference to the accompanying drawings, it is to be understood that the 
invention is not limited to those precise embodiments, and that various other changes and 
modifications may be made by one skilled in the art without departing from the scope or 
spirit of the invention. 
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